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The corrections induced by light quark masses to the current algebra result for the 
7T° lifetime are reexamined. We consider next-to-next-to- leading order corrections and 
we compute all the one-loop and the two-loop diagrams which contribute to the decay 



■ amplitude at this order in the two-flavour chiral expansion. We show that the result 

42 is renormalizable, as Weinberg consistency conditions are satisfied. We find that chiral 

logarithms are present at this order unlike the case at next-to-leading order. The result 
could be used in conjunction with lattice QCD simulations, the feasibility of which was 
recently demonstrated. We discuss the matching between the two-flavour and the three- 
flavour chiral expansions in the anomalous sector at order one-loop and derive the relations 
between the coupling constants. A modified chiral counting is proposed, in which m s 
counts as 0(p). We have updated the various inputs needed and used this to make a 
phenomenological prediction. 
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1 Introduction 

X 

The close agreement between the current algebra prediction for the lifetime of the neutral 
pion and experiment is one of the two compelling experimental signatures, together with 
the Nambu-Goldberger-Treiman relation, for the spontaneous breaking of chiral symmetry 
in QCD. There is an ongoing effort by the PrimEx collaboration [1] to improve significantly 
the accuracy of the lifetime measurement, which is now around 8%, down to the l%-2% 
level. This motivates us to study the corrections to the current algebra prediction. 

Starting with the detailed study by Kitazawa [2], this problem has been addressed 
several times in the literature [3-11]. The approach used in ref. [2] was to extrapolate 
from the soft pion limit to the physical pion mass result using the Pagels-Zepeda [12] 
sum rule method. This was reconsidered in ref. [10] who implemented a more elaborate 
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treatment of 7r° — r\ — rj mixing and also recently in ref. [11]. The latter work shows some 
disagreement concerning the size of the rj meson contribution in the sum rule as compared 
to earlier results. 

In this paper, we reconsider the issue of the corrections to the current algebra result 
to the 7T° — > 27 amplitude from the point of view of a strict expansion as a function 
of the light quark masses. This is most easily implemented by using chiral Lagrangian 
methods (see e.g. [13] for a review). The same framework also allows one to compute 
radiative corrections [14]. We believe that it is somewhat easier to control the size of the 
errors in this kind of approach, which is important for exploiting the forthcoming high 
experimental accuracy. Another interest in deriving a quark mass expansion is the ability 
to perform comparisons with lattice QCD results where quark masses can be varied. The 
feasibility of computing the 7r° to two photon amplitude in lattice QCD has been studied 
very recently [15]. 

A priory, it is expected that one can make use of SU(2) ChPT, i.e. expand as a 
function of m u , rrid without making any assumption concerning m s (except that it is 
heavier than m u , m^). In SU{2) ChPT it is often the case that chiral logarithms provide 
a reasonable order of magnitude for the size of the chiral corrections. This is the case, 
for instance, for the tttt scattering amplitude [16,17]. It was observed in refs. [3,4] that 
there was no chiral logarithm in the next-to- leading order (NLO) correction to the 7r° 
lifetime once the amplitude is expressed in terms of the physical value of F n . We have 
asked ourselves whether chiral logarithms are present in the NNLO corrections. At this 
order, the coefficient of the double chiral logarithm depends only on F n . Depending on its 
numerical coefficient, such a term could modify the NLO results. In order to obtain this 
coefficient it is, in principle, sufficient to compute a set of one-loop graphs containing one 
divergent NLO vertex [18]. For completeness, we will perform a complete calculation of 
the twodoop graphs as well. This is described in sec. El 

In the framework of two-flavour ChPT, however, one faces the practical problem that 
the polynomial terms in m u , at NLO involve a number of low-energy couplings (LEC's) 
which are not known. We will show that it is possible to make estimates for the relevant 
combinations, and then make quantitative predictions for the 7r° decay, under the minimal 
additional assumption that the mass of the strange quark is sufficiently small, justifying 
a chiral expansion in m s . We will obtain the first two terms in the m s expansion of 
the NLO SU(2) LEC's. The result can be implemented in association with a modified 
chiral counting scheme, in which m s is counted as 0(p), which respects the hierarchy 
Tn u , rrifi <C m s . This leads to simpler formulas than previously obtained. Finally, we will 
update all the inputs needed to compute the lifetime. 

2 Leading and next-to-leading orders in the SU (2) expansion 

In the odd-intrinsic-parity sector, the Lagrangian of lowest chiral order has order p 4 , it is 
the Wess-Zumino [19] Lagrangian, C wz , which form is dictated by the ABJ anomaly [20]. 
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Writing the tt° — ► 7(^1)7(^2) decay amplitude in the form 

T = e 2 e(e* 1 ,k 1 ,e* 2 ,k 2 )T , (1) 
a tree level computation of the pion decay amplitude gives the well known result 

where F is the pion decay constant in the two- flavour chiral limit m u = = 0. At 
leading order one can set F = F v in eq. ([2]). According to the Weinberg rules [18] for 
ChPT, the NLO corrections are generated from: 

a) One-loop diagrams with one vertex taken from C wz and other vertices from the 
0(p 2 ) chiral Lagrangian. These diagrams were first computed in refs. [3,4]. 

b) Tree diagrams having one vertex from C wz and one vertex from the 0{p i ) chiral 
Lagrangian. 

c) Tree diagrams from the 0(p®) Lagrangian in the anomalous-parity sector, £^§y 

The classification of a minimal set of independent terms in this Lagrangian was initiated 
in refs. [21,22]. We will use here the result of ref. [23] who further reduced the set to 23 
terms in the case of three flavours and to 13 independent terms in the case of two flavours 
(this result was also obtained in ref. [24]). The list, in the case of two flavours, is recalled 
below: 

jC Sv>=2 = e Q/3 ^{ cf (X+ [f-fa,, U Q Up\) +cl V (X-{U t ,u,U a U f} }) +cfi(X-f+^Uaf3) 
+cfi(X-f-^f-al3) +cfi(x+[f+^J-a/3}) (f+fiw) {x~U a Up) 

+cY i (f+vv)(f+apX-) +cYi{f+^)(f +a p){x-)+Cv iif+j^ihluaUp) 

+ClO^(/+ ,j}(f-^U a Up) +cY 1 (f + ^)(f+ 1 ah' y i3 ) +C^(/ + ^)(/7 a /_ 7/3 ) 
7/Lt)(/W"/3}} • (3) 
The relations between the bare and the renormalized couplings may be written as [23] 

JV _ Wrr s , W_M_I__ (A) 

c % -c, W+r, t 167r 2 (d _ 4) W 

with log(c) = — (log(4-7r) — 7 + l)/2 as usual in ChPT (note that the couplings cf T have 
dimension (mass) -2 ). The coefficients vanish for % = 1...5 and the remaining ones 
read [23] 

rj^ = 3a, v Y = 3a, r]f = -|a, r$ = 6a 
< = -18a, < = 12a, r$ = 0, rj% = -12a , 

with 

a = 1/(384^ 2 F 2 ) . (6) 



3 



The above results for 'qf were obtained by using, in the ordinary sector at p 4 , the chiral 
Lagrangian term proportional to I4 which differs from the form originally used in ref. [25] 

£ (orig) = ^ {u ^_ } (7) 
by a term proportional to the equation of motion 

A 4 = 4 oris) + ^(x-(v^-lx-)). (8) 

If one uses >c|° rig ^ then, in the odd-intrinsic-parity sector, the coefficients with labels 6,7 
and 8 are modified to cf^ [7,26]. The relations between and are easily worked out 
by performing a field redefinition, 

~w _ w N c I4 



c 6 — c 6 



128vr 2 F 2 



07 — C7 T" 



256vr 2 F 2 

d w = c w _JY^_il . (9) 
8 8 512vr 2 F 2 v ' 

In the present work, we use the original £;°" g ^ in our calculations but we will express 
the final result in terms of c^ r rather than c^ r , making use of the relations ([9]) (which 
will prove slightly more convenient below when we perform a matching with the SU (3) 
expansion). 

Returning to the tt° decay amplitude, the contributions from the one-loop Feynman 
diagrams can be shown to be absorbed into the re-expression of F into F v [3,4], the 
physical pion decay constant at order p , such that the decay amplitude including the 
NLO corrections reads 

rr 1 J 1 t 16 2 / a Wr a Wr , Wr\ 



F n I 4vr 2 3 
64 



g B(m d - m u )(5c^ r + c™ r + 2c 8 ^) j , (10) 

where B = — lim mii=md= o (uu)/F 2 and m 2 denotes the mass squared of the neutral pion 
which, at this order, is equal to M 2 = B{rrid + m u ) . Eq. (I10p shows that the decay 
amplitude at NLO receives a contribution proportional to the isospin breaking mass dif- 
ference rrid — m u . As can be seen from eqs. (JSJ) the two combinations of chiral couplings 
which enter into the expression of T/vlo are finite. The expression of Tlo+nlo there- 
fore involves no chiral logarithm. The chiral corrections to the current algebra result are 
purely polynomial in m u , and are controlled by four coupling constants from eq. ([3]). 
In order to estimate quantitatively the effects of the NLO corrections, we will show below 
that it is useful to express these couplings as an expansion in powers of the strange quark 
mass. Before doing so, let us now investigate the presence of chiral logarithms, which 
could possibly be numerically important, in the NNLO corrections. 
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3 7T° decay to NNLO in the two-flavour expansion 

We must calculate now a) the one-loop Feynman diagrams with one vertex involving an 
NLO chiral coupling, either li or and b) the two-loop Feynman diagrams with one 
vertex taken from the LO Wess-Zumino Lagrangian and the other one taken from the 
0(p 2 ) chiral Lagrangian. It is convenient the use the following representation for the 
chiral field 



.T • 7T 

U = a + i — , a 



TT 2 



• " (11) 

p ' V F 2 

(since, in this representation, there is no "fin vertex at LO). At the order considered, all the 
reducible diagrams are generated from wave- function renormalization. The expression for 
the WF renormalization constant Z (corresponding to (jlljl ) was first given by Biirgi [27], 



i 

Z2 



1 



Tm 

2F 2 ~*~ F 4 



\t 2 m + ^(r z + f 2 M Q z - t M hQl 



+ B 2 {m d - m u ) 2 -8F 2 (c 7 + c 9 ) + t M h 



(12) 



with 



dT M 



(4vr)f ' M dM 2 



, d = 4 + 2w . 



(13) 



We have indicated explicitly here the contributions proportional to (m<2 — m u ) 2 for com- 
pleteness because isospin breaking contributions play an important role for the 7r° decay 
amplitude. We will also need the expression for the chiral expansion of F n at order p 6 
(from [28]) 



r F + f M Q F -t M J2UQ 



8B 2 (m d -m u ) 2 
H ^ (C7 + eg J 



F 2 



The numerical parameters Q z and Q F which appear above read 



(14) 



Q z = —(96 -464w + l 185w 2 ), Q F = —(240 - 656u; + 1125w 2 ) (15) 

96 192 



and will need the following relations obeyed by the numerical parameters Qf and Q\ 

1 



—Qi Q2 



T Q2 > Qi -Qi ~ 2 > Qi - 2(1 + w ) 



(16) 



Finally, the entries rz and rp in eqs. (|12fl . (j!4H represent combinations of coupling con- 
stants from the 0(p ) chiral Lagrangian. The 7r° amplitude involves the combination 
rz + 2rp which is expressed in terms of a single p 6 coupling, called cq in the classification 
of ref. [29] 



r z + 2r F = —6&F c 6 



(17) 
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(a) (b) 



(c) (d) 




The two-loop one-particle irreducible diagrams which one must compute (using the rep- 
resentation (jlip ) are shown in fig. [TJ It turns out to be possible to express all of them 
analytically in terms of known special functions by combining the methods exposed in 
ref. [30] with integration by parts methods. We give the results corresponding to the two 
diagrams (f) and (g), which are the most difficult ones, in appendix 1. 

Collecting all the pieces together, we find that the expression for the NNLO contribu- 
tion to the 7T° decay amplitude into two photons has the following expression, 

M 4 / 1 



M 4 
+ 16^ 2 F 4 



^ {2cf r + 4cf r + 2c^ r + 4cY r - cf x r ) 



256?r 4 

+ 32M2 ^ d F ~ mu) [-6cT r ~ Hcf r + 6cf r - 12cf - cY r - 2cY r ) 
I M4 ,V I M2B ( m d- m u) x | B 2 (m d -m u ) 2 x (18) 

where L n represents the chiral logarithm 

9 

TO 

L^ = log^ (19) 
and A + , A_, A can be expressed as follows in terms of renormalized chiral coupling 



6 



constants, 



A + — — 2 



A. 

A_ 



-Mf + 



1 



5127T 4 



QQO A \ " 

-— --C(3) + 3V3Cl 2 (7r/3)J 



16 



+ f^ 2 [^(c3 W/r + cD+^(-4c3 

^ [d^(Ai) + F% {bcf r + cY r + 2cf r )] 
9 

= d^ZQi) - l2%F 2 l 7 {cf r + cf r ) . 



(20) 



Here, the notation (i 1177 * refer to combinations of couplings from the NNLO Lagrangian 
(i.e. of order p 8 ) in the anomalous sector. 

A few remarks are in order concerning this calculation. First, concerning non-local 
divergences, i.e. terms of the form M 4 log(M 2 )/(d — 4), we have verified that those which 
are generated from the two-loop diagrams are cancelled exactly by those generated from 
the one- loop diagrams proportional to Zj, as expected from the Weinberg consistency 



conditions. The divergences that are left are proportional to M , M (md- 



■m„ 



and (m,d- 



m u ) . They are cancelled by the contributions, at tree level, from the chiral Lagrangian of 
order p 8 in the anomalous sector. We have denoted the three independent combinations 
of 0(p s ) chiral couplings by , d_ and d^_. Our calculation shows that the relation 
between these and the corresponding renormalized combinations must be as follows, 



d 



w _ M 2(d - 4) 



F 2 



df» - A 2 (- 
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A 



nil 



ill 



l 



• r 

:1a 



53 



4608vr 2 



+ 16tt 2 F 2 



-4cY r 



2c 



Wr 



4c, 



Wr 



2o 



Wr 



iW 



(c/z) 



2(<2-4) 



d 



F 2 



[d^ r (^) - AF 2 



23cf r + 18cf r 



36c 



Wr 



JVr 



F 2 



7T Z 



(211 



Eq. (|18|) shows that chiral logarithms are indeed present at NNLO. The coefficient of the 
dominant one, as can be shown quite generally, depends only on F. The coefficient of 
the subdominant chiral logarithm has one part depending only on F and another one 
depending on the NLO chiral couplings cf r - From a numerical point of view, the con- 
tribution from the dominant chiral logarithm turns out to be very small, of the order of 
a few per mille. This lack of enhancement could indicate a fast convergence of the chiral 
perturbation series. In this respect, the detailed formula (fT8|) could be used in association 
with results from lattice QCD simulations, in which the quark masses m u , are larger 
than the physical ones and can be varied. This would allow one to determine the relevant 
combinations of chiral couplings. In the following section we discuss an alternative, more 
approximate method, to estimate these combinations. 
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4 Chiral expansion in m s 



From now on, we assume that the mass of the strange quark is sufficiently small, such 
that the chiral expansion in m s is meaningful. One can then calculate the 7r° lifetime 
using the three-flavour chiral expansion. Instead of doing so directly, as it remains true 
that m u , m,i « m s, it is instructive to start from the SU{2) expression, eq. (fTSJ) and 
perform a chiral expansion of the couplings cf^ as a function of m s . A priory, one expects 
expressions of the following form to arise 

cY r = ^ + U + £ ^ C « T + <Uog^^ ) + 0{m s ) , (22 



m. 




where C^ r , a = 1 ... 24 are the coupling constants of the NLO Lagrangian in the anoma- 
lous sector in the SU(3) expansion [23] and Bq = lim ms= o B. Analogous expansions were 
established in ref. [31] for the SU(2) couplings B, F and l\. This problem was reconsidered 
recently in ref. [32] in which the NNLO terms in that expansion have been derived. Also in 
ref. [33] the m s expansions of the SU{2) LEC's in the electromagnetic sector were studied. 
In order to generate such expansions one can work in the SU{2) chiral limit m u = rrid = 0, 
compute sets of correlations functions having SU (2) flavour structure in both the SU{2) 
and the SU (3) chiral expansions and equate the expressions. The authors of ref. [32] have 
shown how to perform this matching at the level of the generating functionals. In the 
SU{3) generating functional, one must use external sources s, p, v^, a M which correspond 
to those used in the SU(2) functional embedded into 3x3 matrices. Since there is no 
source for strangeness, the classical SU(3) chiral field involves the three pions ir a and the 
7] field but no kaons 

U d = exp exp — — (23) 

(Fq being the pion decay constant in the three-flavour chiral limit). Using the equation of 
motion one can express the field rj c i in terms of an SU{2) chiral building-block [32,33] 

// "' i(X-)(-7^-5)+0( P 4 ). (24) 



V3F V 16m SJ B 

The terms proportional to tj c i thus generate contributions proportional to l/(m s B). These 
can be also seen as resulting from 7] meson propagators in tree diagrams. Besides, eq. (|24D 
shows that r] c i counts as 0(p 2 ) in the SU{2) counting. Inserting U c i from eq. (j23[) in the 
SU(3) Wess-Zumino action and expanding to first order in r\ c i we obtain, 

+ ^</+ Q /3> (/+/-> " \i{f-apf-^)\ ■ (25) 

This allows one to deduce the leading terms, which behave as l/m s , in the expansion 
of the couplings c-^ r . Next, the terms proportional to (m s )° are generated from three 
sources. 
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1) Prom the SU(3) Lagrangian C w , by inserting U c i (with rjd set to zero), which gives 
contributions proportional to LEC's Cf^ T '. 

2) From one-loop irreducible graphs with one vertex taken from the Wess-Zumino action 
and having one kaon or one eta running in the loop. 

3) From corrections to the r/ pole contributions stemming from tadpoles or from vertices 
proportional to the 0(p 4 ) couplings L». 



The results are presented in eqs. (f2"6"|) below and (|A-4|) in the appendix. 

Let us now examine the applications of this exercise to the problem of the 7r° lifetime. 
As seen in sec. [2] the NLO corrections involve two independent pieces, one proportional 
to mi and one to B(rrid — m u ), and they are controlled by two combinations of the four 



couplings dj r , Cy r , c§ r and C\\ ■ For these, we take into account the first two terms in 
the m s expansion which read 



Wr 
c 3 



Wr 
c 7 



Wr 
c 8 

Wr 
c ll 



c + Cf r + 3C^ r + 0(m 



Wr 



Wr 



2 

3 

2 C0 



~C% r + 0(m s 



(26) 



where 



and 



co 



32tt 2 



1 



l6Bm, F, 



L K = log ■ 



o 



m,B n 



+ -^(3L r 7 + Ll 



1 IT 2 r 

512^ ( ^ + 3^ 



L„ = L K + log ■ 



(27) 



(28) 



fi" ' 3 

At this point, one observes that by using the m s expansion, we have expressed four SU{2) 
couplings in terms of three SU{3) ones. This might look as a modest improvement. 
Fortunately, the combinations relevant for the tt° lifetime at NLO actually involve only 
two couplings C^ r , C^ r while C^ r drops out. 



Let us now consider the terms proportional to ra\ and m 2 (m^ 



One can see 



from eq. (|18p that they involve four more LEC's, c^ 7 ", c^ 1 ", c^ r , c§ T ■ It makes sense here 
to retain only the part of these LEC's which are dominant in the m s expansion, i.e. the 
l/m s part, 



Wr 



CO, 



Wr 
C4 ~ 



1 



:C0, 



cf - 0, 



Wr 

c 6 - 



1 



-co, c 



512^5^, 



(29) 



and we perform a similar approximation in eq. (I26p . We will also retain the part involving 
the LEC Co^ as it will appear that the size of this coupling is comparable to that of the 
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l/m s terms. Inserting the m s expansions (I26p and (|29p . in the SU(2) chiral expansion 
of the 7T° decay amplitude (fTHj) we obtain the following expression 



T, 



4 



1 1 64 2 rWr , I m d -m u 



16vr 2 



+ 32B(m d - m u ) 



3 G 7 + 4C 8 ^1 3^ 



3 m 2 

2 16^ 2 F 2 * 



(30) 



1 / 12 



1 



ml 



24tt 2 Vl67r 2 F2 L7r 



4.1 A modified SU(3) chiral counting 

Some comments are in order concerning eq. (|30[) . In particular, one expects it to be related 
to the formula that one can compute starting from SU(3) ChPT. Such a computation was 
performed, e.g. in ref. [7]. In SU(3) ChPT m u , m d and m s are counted on the same 
footing, 

f^ui iTid ~ m s ~ 0(p 2 ) [standard SU(3)] . (31) 

In the physical situation, however, m u , « m s . For processes which involve only 
pions this can be accounted for by adopting the following modified counting, 

m u , m d ~ 0(p 2 ), m s ~0(p) [modified SU(3)] . (32) 

The formula (|30l) for the 7r° lifetime can be argued to be a consistent expansion in this 
modified counting. One notes first that all the corrections must be proportional to m u , m d 
since the starting point is exact in the SU(2) chiral limit. The formula (|30p includes the 
leading corrections of order p (which must be proportional to m u /m s , m d /m s ) as well as 
the subleading corrections of order p 2 (which must be proportional to m u , m d ). It also 
includes the corrections of order p 3 which are logarithmically enhanced (which must be 
proportional to m u m s , m d m s multiplied by log(m u -\-m d ) as well as the corrections of order 
p which are double logarithmically enhanced. Obviously, by retaining logarithmically 
enhanced terms at a given order instead of the full set of terms, one introduces a chiral 
scale dependence into the amplitude. Clearly, one should use a value of the scale of the 
order of the kaon or the eta mass for this approximation to make sense. Finally, we have 
verified that, starting from the expression for the amplitude in standard SU (3) at NLO 
obtained in ref. [7], and expanding in powers of m u /m s , m d /m s one recovers exactly the 
terms of order p, p 2 and p 3 log(p 2 ) in the modified SU(3) expansion (|32p . In practice, 
the expression (|30|) is somewhat simpler than the standard SU(3) NLO expression and 
contains the double logarithm term. The latter turns out to be numerically small so that 
the two expressions are essentially equivalent in practice. In order to derive a numerical 
prediction from eq. (j30|) one needs inputs for: F w , (rrid — m u )/m s , B(m d — m u ) and , 
Cg^. We will give an update on the determination of these quantities in sec. [5] 

In addition to the chiral corrections induced by the quark masses, one should also take 
electromagnetic corrections into account. These have been considered in ref. [7], where 
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the correction terms of order e 2 and of order e 2 (m u + md)/m s have been computed. Here, 
it is consistent to retain only the term of order e 2 , its expression in terms of Urech's chiral 
couplings [14] is recalled, 

= + Kl) + 2KI - Kl - ^ + Jtf) + ^^(5 + 4L W + L K )} . 

(33) 

This term is defined such that the 7r° — ► 27 amplitude is expressed in terms of F n o which 
is the neutral pion decay constant in pure QCD and m 2 which is the physical neutral pion 
mass (i.e. including EM corrections). 



5 Phenomenological updates 

Let us now update the various inputs needed to calculate the numerical prediction for the 
7T° lifetime in ChPT. 



1) F v 



An obviously essential input here is F n , the value of the pion decay constant. 
Marciano and Sirlin [34] have evaluated the radiative corrections in the process 
7r + — > decay rate such that it is expressed in terms of F w + the charged pion 

decay constant in pure QCD. In pure QCD the difference between F^o and F n + is 
quadratic in the quark mass difference — m u and can be expressed as follows in 
ChPT, 



F„ 



F„o 



1 



QCD 



B 2 (m d - m u ) 



-16c&(/i) 



h 



16vr 2 



1 + log 



ml 



~ 0.7 10" 4 . 



A rough numerical evaluation has been made by using leading order l/m s estimates 

P2 q / 7^2 n 2 



h 



F~ 



8Bm< 



^9 



2 \Bm< 



(34) 
lates 

(35) 



Eq. (|34p shows that the difference between F n + and F w o is negligibly small for our 
purposes, and we will ignore it. In the expression of ref. [34] for the radiative 
corrections, one constant term, called C\, was left undetermined. Matching with the 
ChPT expansion of the ir + decay rate at 0(p 4 ) one can express C\ in terms of chiral 
logarithms and a set of chiral couplings [35]. The latter can then be estimated using 
chiral sum rules and resonance saturation [36]. Using these results and the updated 
value of Vud from ref. [37] 

V ud = 0.97418(26) , (36) 



we find 



F n = 92.22 ± 0.07 MeV . 



(37) 
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B{m d - mu), (m d - m u )/m s , 3L 7 + L r 8 : 

Because of the Kaplan-Manohar invariance [38] it is not possible to determine in- 
dependently the quark mass ratios and the couplings L7, L% in ChPT using low- 
energy data. One may use an input from lattice QCD, e.g. on the quark mass ratio 
r = 2m s /(m u + m^). Using the results obtained in ref. [39] as well as those from 
other recent QCD simulations which are collected in table XVI of that reference and 
averaging, one can deduce 

2/ "- 28.0 ±1.5. (38) 



m u + m d 



Using this input for r, we may treat terms linear in the quark masses in NLO ChPT 
expressions as follows, 

(m u + m d )B ~ ml, m s B = . (39) 

The value of the LEC combination 3Lj + Lg, can be deduced using r and standard 
0(p i ) ChPT formulas for the pseudo-scalar meson masses [31] 

3L 7 + L r 8 (n) = (0.10 ± 0.06) 10~ 3 (fx = M v ) . (40) 

Concerning the quark mass difference m d — m u , we will use the recent determination 
made in ref. [40]. It is based on the 77 — ^ 37r decay amplitude which is an isospin 
breaking observable with very small sensitivity to electromagnetic effects [41,42]. 
The amplitude has been computed at order p 6 in ChPT by the authors of ref. [40] 
and they deduce the following result^ , 

R =. = 42.2 (41) 

m d - m u 

(with m = (m u +m d )/2). No figure for the uncertainty is given. We have estimated it 
by noting that the main source of uncertainty in this result comes from the unknown 
values of the coupling constants CJ from the 0(p e ) Lagrangian. For these couplings, 
it was shown that simple resonance models are sometimes misleading [44] because 
of their strong scale dependence. We have estimated the order of magnitude of 
the uncertainty by taking the difference between the value of R obtained from a p 6 
calculation and the value obtained from a p 4 calculation and dividing by two, which 
gives 

AR ~ 5 . (42) 



1 An alternative evaluation of R can be made based on the K + — K° mass difference. As one can see 
from table 6 of ref. [40] this method tends to give values of R smaller than eq. (|4ip . The calculation of the 
K + — K° mass difference in ChPT, however, has uncertainties related to the couplings d and also from 
estimates of the electromagnetic contributions, beyond the Dashen low-energy theorem, which have some 
model dependence. One could also use isospin violation in K13 form factors. For an updated discussion of 
these effects see [43]. 
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Using dHJ and (|I2]), we obtairH 



— — — = (2.29 ± 0.23) 10" 2 , B(m d - m u ) = (0.32 ± 0.03) M\ . (43) 
m s 

3) C^: 

This constant obeys a sum rule in terms of the form factor associated with the 
photon-photon matrix element of the pseudoscalar current [2]. A simple resonance 
saturation approximation in this sum rule gives a relation between and the 
7r(1300) mass and its couplings to the pseudoscalar current (d m ) and to two photons 
GfcrO [7] 

- ^ • (44) 

7r' 

Recent experimental data by the Belle collaboration has confirmed the extreme 
smallness of the coupling of the 7r(1300) meson to two photons [46] 



+ 2o 



< 72 eV . (45) 



The validity of the resonance saturation approximation in this case might be ques- 
tioned since, in the sum rule, , could pick up more important contributions from 
energies higher than the mass of the 7r(1300) resonance. There has been several 
attempts at estimating this high energy contribution to in the literature: Us- 
ing a quark-hadron duality picture, Kitazawa [2] argue that this contribution arises 
from a triangle diagram and should thus be proportional to the constituent quark 
mass (this result was applied to r/ decay in ref. [47]). In QCD, one expects the 
constituent quark mass to be momentum dependent (see e.g. [48]) and to decrease 
at high momenta, which is not taken into account in this evaluation. A calcula- 
tion of the triangle diagram in the NJL model was performed in ref. [49]. As this 
model implements a momentum cutoff, however, it rather concerns the low-energy 
rather than the high-energy contribution to Cy . An alternative idea was proposed 
in ref. [50] based on a minimal resonance saturation modelling of the three-point 
function VVP and enforcing a correct asymptotic matching to the OPE expansion 
of this three-point function. The result, unfortunately, cannot be shown to remain 
stable under inclusion of more resonances. None of the estimates, finally, appear to 
be quantitatively very compelling. It seems however quite safe to assume that the 
coupling CY should be suppressed, say by one order of magnitude, as compared to 
the coupling . Indeed, in an analogous sum rule representation, Cg 17 picks up a 
strong contribution from the r/ resonance. We will therefore take 

\Cf | < 0.1 |Cf | . (46) 



2 In ref. [45] a determination of the quantity Bo(md — m») from 77 — > 3tv° was proposed, based on using 
both the decay rate and the slope parameter a, obtaining Bo (mj — m u ) = (0.25 ± 0.02) M 2 • This appears 
somewhat smaller than the result in eq. (|43[1 but one should keep in mind that the ratio Bo/B, while 
expected to be close to one, is not accurately known. 
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4) Cf : 



Having assumed the validity of SU(3) ChPT, together with the result (|46[) of the 
above discussion on , one can determine from the experimental information 
on the r\ — > 27 decay width. According to the PDgH [53] 



77^27 



0.510 ± 0.026 keV , 



(47) 



while the corresponding amplitude computed in ChPT, including LO and NLO con- 
tributions, reads 



F 64 
i ^r(l + x,)- T m^ 



w 



256, 

H (r 

3 V 



l)ml^-CY + (' 



w 



+ 0(m 



where x v encodes isospin breaking effects 



Vs(-ei + e 2 (5 v - 6x)) ~ -0.023 , 



(48) 



(49) 



using notations and results from refs. [31] and [7]. We need an input for in 
eq. ([4"8"]h Up to corrections quadratic in m s , F n is linearly related to F n and Fx [31], 



AF K - F v 



+ 



rat 



967r 2 F w 



2(r + 1) log 



2(2r + 1) 
3(r + 1) 



log 



2r + 1 



+ 0(m^ 



(50) 



The review in ref. [54] quotes the following result for Fk from averaging over recent 
experiments on 7T;2 and K[ 2 decays 



F K Vu 



0.27599(59) . 



F n V ud 

Assuming exact CKM unitarity we can deduce Fk and then F„ 



Fk = 109.84 ± 0.63, F v = 118.4 ±8.0 



(MeV) . 



(51) 



(52) 



The error on F„ is dominated by the 0(m 2 s ) contributions in eq. (I50p . We have 
estimated that it should be smaller than the 0(m s ) contribution by a factor of 
three. Finally, using these results in conjunction with eqs. (|47|) (|48|) we determine 
the coupling 

Cf = (0.58 ± 0.20) 10~ 3 (GeV~ 2 ) . (53) 

We have estimated that the uncertainty stemming from unknown 0(m? s ) chiral cor- 
rections in the rj decay amplitude to be of order 30% compared to the 0(m s ) cor- 
rections. 



The numerical results for the current algebra amplitude and the corrections according 
to the modified chiral SU(3) counting, using the updated inputs presented above, are 

3 The PDG now rejects the Primakoff experiment [51] which gave a smaller result. A re-discussion of 
that experiment has recently appeared [52]. 
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CA 


0(p) 


0(p 2 ) 


0(e 2 ) 


0(p 3 log p) 


0(p 4 log 2 p) 


7.76 


0.09 


0.29 


-0.05 


0.005 


-0.004 



Table 1: Current algebra contribution to the 7r° — > 2j decay width (in eV) and corrections 
of various chiral orders using the modified SU (3) counting. 



collected in table [TJ One remarks that the 0(p 2 ) contribution is larger than the 0(p) one. 
This is induced by the size of the LEC Cf. Expressed as a sum rule, is dominated 
by the r/ contribution, which can be written [7] 

Cf ~ , (54) 

where Mi is the mass of the if in the chiral limit. In the large N c limit one has, 

" 128^' dm " 2^6' 8 ~ 256vr2 M 2 ,' (&5) 

The enhancement of can then be understood, qualitatively, as a large N c effect. In 
practice, the value of that one can estimate using the resonance saturation formula (|54|) 
agrees reasonably well with the one deduced from a ChPT expansion of the rj — > 2j 
amplitude^ (eq. (|53p ). The enhancement of the 0(p 2 ) contribution is therefore a well 
understood effect and does not signal a breakdown of the expansion. Table [T] shows that 
the logarithmically enhanced contributions of order p 3 log(p) and p 4 log 2 (p) are quite small 
in practice and tend to cancel each other. Finally, the prediction for the tt° decay width 
reads, 

r 7r o_ >27 = (8.09 ± 0.11) eV . (56) 

The two main sources for the uncertainty are: — m u (±0.05) and Cg 17 (±0.098). We 
have added the errors in quadrature. Compared to ref. [7] the main modification in the 
input is the value of the rj — > 2~f width in the PDG. The branching fraction for the 27 
decay mode is (98.798 ± 0.032)% [53] (the most sizable other decay being the Dalitz mode 
7T° — ► 7e + e~, for review see e.g. [55]). Our result, eq. (|56|) . then corresponds to the 
following value for the 7r° lifetime 

r^o = (8.04 ±0.11) 10~ 17 s . (57) 



6 Summary 

In this paper, we have reconsidered the chiral expansion of the 7r° — ► 27 amplitude. At 
first, we have focused on the two-flavour expansion. We have considered the expansion 
beyond the known NLO (which we have expressed in terms of the coupling constants 

4 Our result disagrees with ref. [11] in which the corresponding contribution is smaller by one order of 
magnitude. 
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introduced in ref. [23]). We have computed all the loop graphs which contribute at NNLO. 
As expected, we found that the divergences are renormalizable by Lagrangian terms of 
chiral order p s in the anomalous sector. We found that chiral logarithms are present at 
this order. For physical values of the quark masses m u , these NNLO corrections turn 
out to be negligible. Even the terms enhanced by logarithms are numerically very small 
in practice. Our final expression (eq. (fTBj) ) could be useful in association with lattice 
QCD simulations in which unphysical quark masses can be used. This would provide a 
direct evaluation of the SU(2) couplings. As an interesting application, one could deduce 
(using also experimental data such as from PrimEx) a precision determination of F n 
uncorrelated with the value of V u d- Such simulations have not yet been performed for 
correlation functions in the anomalous sector, but this would be of obvious interest. 

In order to perform a more detailed phenomenological analysis at present, it is possible 
to enlarge the chiral expansion from SU(2) to SU(3). This allows one to derive some 
information on the SU (2) coupling constants. We have derived the expansion of the SU (2) 
couplings c^ r as a function of m s up to 0{m s ) and inserted this result into the SU{2) 
expansion formula. The leading, l/m s terms in this expansion, reflect the influence of tt°— r\ 
mixing. We then implemented a modified chiral counting in which m s is counted as 0(p) 
rather than 0(p 2 ). This counting accommodates the fact that m u , are significantly 
smaller than m s . The formulas obtained in this way are somewhat simpler and easier to 
interpret than those obtained in the usual chiral counting but the numerical results are 
essentially identical. 

We have updated the inputs to be used in the chiral formula. A key input is the 
value of F n , the pion decay constant in pure QCD. Another important input is the value 
of the 7] — > 2j decay width, which we use to determine the value of the SU(3) LEC 
Cg^. In the chiral approach, this LEC encodes the effect of t] — rf mixing. Our result 
agrees well with that of approaches which account for r/ — r( mixing explicitly, using 
large N c arguments in addition to chiral counting [2,9, 10]. The overall uncertainty is 
dominated by the unknown terms of order p 3 , i.e. proportional to m u m s , m^m^ in the 
chiral expansion. As a final remark, we note that F n is determined from the weak decay 
of the 7r + assuming the validity of the standard model. Some recently proposed Higgsless 
variants can accommodate deviations from the standard V — A coupling of quarks to the 
W as large as a few percent [56]. Precision measurements of the 7r° lifetime can provide 
constraints on such models. 
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Appendix I 



We give here the result of our computation of diagrams (/) and (g) in fig. [TJ 



M 4 



11 



F « T U)-^i\ ~ -J 



A 2 + A[L n 



31 \ 1 



31 



1056vr 2 ) + 2 Ln 1056vr 2 Ln + 6144vr 2 



467 1 

" 98304tt 4 J 



(A-l) 



F * T (g) ~ n 2p4 | 



M 4 [ 7 
3 



A + A L 



59 



59 r 1 



1 



5127T 4 



fc/5Cb<f ) - i<(3) 



77 1792^ 2 / 2 * 1792vr 2 w 6144vr 2 
1135 



with 



A 



16^ 2 (d-4) 



(A-2) 



(A-3) 



Appendix II 

We collect below the expansions of the SU (2) couplings c^ r as a function of m s up to 
0(m s ). The notations Lx, L v and cq having been introduced in eqs. (|27D and (I28h these 
expansions read 



l„u/„ 1 1_ 

o 



1 ^,M/»- 3 



c 2 = c + G 4 - -G 5 + -G 6 

Q 

C 3 =--c + G 7 +3C 8 

c 4 =-2 C o + ^9 +^10 



Wr _ „ , n Wr 3 Wr 1 Wr 1 n Wr 

c 6 — — Co -I- o 5 — 2 6 _ 2 _ 2 

c 7 — 2 ~~ 8 4 22 

8 ~~ 4 2 ~~ 8 2 



c 9 — ~~ °13 + °14 "t" 15 ~~ o /on_5\') 7^2 ' 1 



3 1 



2 (32vr 2 ) 2 F 2 

Q 1 

Ifr r~<Wr rvfirr i^iWr /~iWr , / j , -i \ 

-10 — °19 ~~ U 20 ~~ U 21 ~ °22 + 2 (327T 2 ) 2 F 2 V 
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Wr 
c ll 

Wr 
c 12 



Wr 
c 13 



fiWr 
< - y 22 



2C%r + 



(32tt 2 ) 2 F 2 



(L k + 1) 



(A-4) 
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